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Abstract 

We present measurements of time-dependent CP asymmetries in B° — > ooKg, fo(980)Kg, Kgir 
and K + K~Kg based on a sample of 535 xlO 6 BB pairs collected at the T(4S I ) resonance with 
the Belle detector at the KEKB energy-asymmetric e + e~ collider. One neutral B meson is fully 
reconstructed in one of the specified decay channels, and the flavor of the accompanying B meson 
is identified from its decay products. CP- violation parameters for each of the decay modes are 
obtained from the asymmetries in the distributions of the proper-time intervals between the two 
B decays. 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 
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The Standard Model (SM) decribes CP violation in B° meson decays using the complex 
phase of the 3x3 Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix CP asymmetries 
in neutral B meson decays into CP eigenstates / exhibit a time- dependent behavior 

A(At) = S f sin(Am d At) + A f cos(Am d At) (1) 

where Sf and Af are the CP violation parameters, Am^ the mass difference between the 
two B° mass eigenstates, At the difference between the decay time of the signal B° (B°) and 
of the opposite-side B° (B°). The SM predicts that for most of the decays that proceed via 
the quark transitions b — > sqq (q = u, d, s) the relations Sf = — sin 20! and Af — 0, where 
£f = +1(— 1) corresponds to CP-even (-odd) final states, hold to a good approximation 0. 
With physics beyond the SM, these decays may receive significant contributions that depend 
on a phase that is different from the SM prediction. A comparison of the effective sin20i 
values, sin 20i , with sin20i obtained from the decays governed by the b — > ccs transition 
is thus an important test of the SM. 

Among the final states studied, luK s and K s 7i° are CP-odd modes, fo(980)K s is a CP- 
even mode, while K + K~K S is a mixture of both £f = —1 and +1. The SM expectation 
for this latter mode is Sf = — (2/ + — 1) sin 20!, where / + is the CP-even fraction. A 
measurement of / + was obtained using isospin relation |(| with a 357 fb^ 1 data sample and 
gives /+ = 0.93 ± 0.09(stat) ± 0.05(syst). 

Recently, it was found that the direct CP asymmetries in B° — > K + ir~ and B + —> K + tt° 
differ significantly 3] while they were naively expected to be same jiJ]. Using the B° — > 
ii'sTr result, an additional test to understand the situation can be made by comparing the 
measured Af value and the value predicted by a sum rule j^] using asymmetry measurements 
from the other B — » Kir decays. 

Previous measurements of CP asymmetries in b — » sqq transitions have been reported 
by Belle and BaBar 0]. Belle's previously published results of CP in B° — > wKg, 
fo(980)Kg, Kg7t° and K + K~Kg were based on a 253 fb _1 data sample corresponding to 
275 x 10 6 BB pairs. In this report, we describe improved measurements incorporating an 
additional 239 fb" 1 data sample for a total of 492 fb" 1 (535 x 10 6 BB pairs). 

At the KEKB energy-asymmetric e + e _ (3.5 on 8.0 GeV) collider, the T(4S) is produced 
with a Lorentz boost of f}^ = 0.425 nearly along the electron beamline (z). Since the B° 
and B° are approximatively at rest in the T(45) center-of-mass system (cms), At can be 
determined from the displacement in z between the two decay vertices: At = Az/(f3jc). 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K\ mesons and to identify muons 
(KLM). The detector is described in detail elsewhere (Iol |. Two different inner detector 
configurations were used. For the first sample of 152 xlO 6 BB pairs, a 2.0 cm radius 
beampipe and a 3-layer silicon vertex detector were used; for the latter 383 xlO 6 BB pairs, 
a 1.5 cm radius beampipe, a 4-layer silicon detector and a small-cell inner drift chamber 
were used [HI]. 

We reconstruct the following B° decay modes to measure CP asymmetries: B° — > ooK s , 
fo(980)K s , K s n° and K + K-R s . We exclude K + K~ pairs that are consistent with a -> 
K + K~ decay from the K + K~K S sample. The intermediate meson states are reconstructed 
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from the following decays: 7r° — > 77, K s — > 7r + 7r~, uo — > 7r + 7r~7r° and /o(980) — > 7r + 7r~. 
Charged tracks reconstructed with the CDC, except for tracks from Ifg — > 7r + 7r~ decays, 
are required to originate from the interaction point (IP). We distinguish charged kaons from 
pions based on a kaon (pion) likelihood £k(k) derived from the TOF, ACC and dE/dx 
measurements in the CDC. Photons are identified as isolated ECL clusters that are not 
matched to any charged track. 

We identify B meson decays using the energy difference AE = E™ s — E^ m and the 

beam-energy constrained mass M bc = J (^beam) 2 — (pl" s ) 2 , where E^ m is the beam energy 
in the cms, and E c ™ and p c ™ are the cms energy and momentum of the reconstructed 
B candidate, respectively. The signal candidates are selected by requiring 5.27 GeV/c 2 < 
Mb c < 5.29 GeV/c 2 and a mode-dependent AE window. The dominant background for 
the b — > sqq signal comes from continuum events (e + e~ — > uu,dd, ss,cc). We discriminate 
against this using event topology: continuum events tend to be jet-like in the cms, while 
e + e~ — > BB events tend to be spherical. To quantify event topology, we calculate modified 
Fox- Wolfram moments and combine them into a Fisher discriminant [12|. We calculate 
a probability density function (PDF) for this discriminant and multiply it by a PDF for 
cos#b, where 63 is the polar angle in the cms between the B direction and the beam axis. 
The PDFs for signal and continuum are obtained from Monte Carlo (MC) simulation and 
a data sideband, respectively. These PDFs are then used to calculate a signal [background] 
likelihood £ S i g [bk g ], and we impose loose mode-dependent requirements on the likelihood 
ratio 7l s /b = £si g /(£si g + £bk g )- Figures UJ (a)- (1) show the reconstructed variables Mbc 
and lZ s /b after flavor tagging and vertex reconstruction (before vertex reconstruction for the 
decay B° — > K s tt°); the corresponding signal yields are summarized in Table |l] 

TABLE I: Estimated signal yields iV S i g in the signal region for each mode. 



Mode 








LUKg 




-1 


118 ±18 






+1 


377 ± 25 






-1 


515 ± 32 


K+K- 


■K s 


+0.86 ±0.18 ±0.09 


840 ± 34 



We determine Sf and Af for each mode by performing an unbinned maximum-likelihood 
fit to the observed At distribution. The decay rate is given by 

e -\At\/r B o , .J 

V{At) = — <^ 1 + q ■ S f sm(Am d At) + At cos(Am d At) } (2) 

4t b o I l J J 

where r^o is the B° lifetime and the 6-flavor charge q = ±1( — 1) when the tagging B 
meson is a B° (B°). The 6-flavor of the accompanying B meson is identified by a tagging 
algorithm ^| that categorizes charged leptons, kaons and A's found in the event. The 
algorithm returns two parameters: the 6-flavor charge q and r, which indicates the tag 
quality as determined from MC simulation and varies from r = for no flavor discrimination 
to r = 1 for unambiguous flavor assignment. If r < 0.1, we set the wrong tag fraction to 
0.5, and therefore the accompanying B meson provides no tagging information in this case. 
Events with r > 0.1 are sorted into six intervals. 
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AE (GeV) N/l c (GeV/c 2 ) Likelihood Ratio 



FIG. 1: AE, M hc and K s/h distributions for (a, b, c) B° -» uK%, (d, e, f) B° -> / if§, (g, h, i) 
£?° — > -f^57T° and (j, k, 1) 5° — » 2f + lif~lf<y. The solid curves show the fits to signal plus background 
distributions, and the dashed curves show the background contributions. To enhance the signal in 
the AE and M^ c projections, an additional cut on TZ S /^ was applied (> 0.5). 



To the PDF expected for the signal distribution (Eq. |2J), the effect of incorrect flavor 
assignment is incorporated and then convolved with a resolution function i? sig (At) to take 
into account the finite vertex resolution. The wrong tag fractions for the six r intervals, wi 
(/ = 1,6), and differences between B° and B° decays, Awi, as well as the resolution pa- 
rameters are determined using a high-statistics control sample of semileptonic and hadronic 
b — > c decays. 
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We determine the following likelihood for each event: 



(1 " U) J 



/ sig P sig (Ai / )i? S ig(A^ - At') 



+ (1 - /si g )nkg(At , )^bkg(At i - At') 



d(At') 



(3) 



The signal probability / S j g depends on the r region and is calculated on an event-by-event 
basis as a function of M bc , AE and TZ s/h (and M(n + n-n°) for B° -> wifg). For 5° 
fo(980)K%, this fit yields the number of B° -> tt+tt 



candidates that have 7T 



7r~ invariant 

mass within the /o(980) resonance region, which includes other contributions (e.g. B° — > 
p°Kg, K*7r and non-resonant three-body decays) which peak like the signal in AE and 
Mb c distributions. To estimate these peaking backgrounds, we perform a fit to the 7r + 7r~ 
invariant mass distribution for the events inside the AE-M^ C signal region (Fig[2J). 



_ 200 




.2 0.4 0.6 



0.8 1 1.2 1.4 

M(tcV) (Gev/c 2 ) 



FIG. 2: 7r + 7r~ mass distribution for the foK$ events in the AE-M^ signal box (shown here after 
background subtraction). The histogram is the result of the fit whereas the different contributions 
are shown (continuous line for /o(980), dashed for p° and dotted for fx)- 



The PDF for background events, "Pbkg(At), is modeled as a sum of exponential and prompt 
components and is convolved with a sum of two Gaussians, i?bk g - Parameters in 'Pbkg(At) 
and i?bkg for background are determined from a fit to the At distribution for events in the 
AE-M^c data sideband. P Q i(At) is a broad Gaussian function that represents an outlier 
component with a small fraction / Q i. The only free parameters in the final fits are Sf and 
Af, which are determined by maximizing the likelihood function L = YliPi{Ati]Sf,Af) 
where the product is over all events. 

Table HTl summarizes the fit results of sin 20^ and Af. For the B° — > K + K~Kg decay, 
the SM prediction is given by Sf = —(2/+ — 1) sin20'{ ff . The effective sin20i value for this 
mode is found to be +0.68 ± 0.15 ± 0.03±g;^. The third error is an additional systematic 
error arising from the uncertainty of the — +1 fraction. We define the raw asymmetry 
in each At bin by (N q=+1 — N q= _i) / (N q=+ i + N q= _i), where A* r g=+1 (_ 1 ) is the number of 
observed candidates with q = +1( — 1). Figure El shows this asymmetry for good tag quality 
(r > 0.5) events in each mode. The dominant sources of systematic error for Sf in b — > sqq 
modes are the uncertainties in the vertex reconstruction (0.01), in the background fraction 
(from 0.01 for Ks'K to 0.04 in ojK^) and in the background At distribution (0.04 in Ks'K 
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TABLE II: Results of the fits to the At distributions. The first error is statistical and the second 
error is systematic. The third error for sin20| of K + K~ K s is an additional systematic error 
arising from the uncertainty of the = +1 fraction. 



Mode 


sin 2(f>f 


At 


ujK s 


+0.11 ±0.46 ±0.07 


-0.09 ± 0.29 ± 0.06 


foK° s 


±0.18 ±0.23 ±0.11 


-0.15 ±0.15 ±0.07 




±0.33 ±0.35 ±0.08 


-0.05 ±0.14 ±0.05 


K + K~K° S 


±0.68 ± 0.15 ± 0.03ig;i3 


-0.09 ±0.10 ±0.05 



and 0.01 or less in others), and in the resolution function (0.05 for uKs and Ksir ). The 
dominant sources for Af are the effects of tag-side interference 0] (0.04), the uncertainties 
in the vertex reconstruction (0.02), in the background fraction (0.03 for foKg and ooKg and 
< 0.02 for others). For the JqK s mode, additional systematics were included: uncertainties 
from the M(irn) fit (0.06 for Sf) and from the assumption on the CP content of the peaking 
background (0.08 for Sf and 0.04 for Af). For the K s ir° mode, the uncertainty on the 
rare B component is taken into account (0.04 for Sf and 0.02 for Af). Other contributions 
come from uncertainties in wrong tag fractions, lifetime and mixing. A possible fit bias is 
examined by fitting a large number of MC events. We add each contribution in quadrature 
to obtain the total systematic uncertainty. 

In summary, we have performed improved measurements of CP- violation parameters 
sin 20^ and A f for B° -> uK° s , f (980)K° s , K s n° and K + K~K° S using 535 xlO 6 BB events. 
Comparing the results for each individual b — > s mode with those from the B° — > J/ipK° 
decay, we have not observed a significant deviation with the present statistics. 
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FIG. 3: Asymmetries of good-tagged events (r > 0.5) for (a) B° -► wKg, (b) £° -> f (980)K° s , 
(c) i?° — > i^^Tr and (d) B° — > K + K~ Kg. The solid curves show the results of the unbinned 
maximum-likelihood fits. The dashed curves show the SM expectation with the measurement of 
CP-violation parameters for the B° — > J/ipK° mode (sin2c/>i = +0.642 and Af = 0) 
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